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The recent observation of a insulator to metal transition (IMT) _I] in pure LaMnOa at 32 GPa and 
room temperature, well above the Neel temperature (145 K) and below the Jahn- Teller transition 
temperature (780 K), opens the way to a study of the role of the orbital degrees of freedom on the 
electronic structure in a stoichiometric material. 

In this paper we focus our attention in the orbital aspects of the insulator to metal transition. 
We use a model Hamiltonian for the Cg orbitals of Mn that includes the on site Coulomb repulsion 
U , the hopping t, and its dependence with pressure. In order to include in an appropriate way the 
strong correlations induced by the dominant electron-electron interactions we introduce auxiliary 
fields (Slave Bosons, SB) to the description of the low energy states. We use the O-Mn distance (d) 
dependence of t according to [H and the pressure-d relation from the experimental data to describe 
the evolution of the electronic structure with pressure. 

Our results confirm and make transparent the conclusion reached in previous ab-initio calculations: 
the inclusion of the Coulomb energy is necessary and constitutes an important factor enhancing the 
orbital polarization in these compounds. 

I. INTRODUCTION 

The discovery of colossal magneto-resistance in manganese perovskytes, of relevance to spintronics and other techno- 
logical applications, has opened a new field to the study of highly correlated systems (HCS). In the three dimensional 
compounds like Lai_a;Sra;Mn03 or in the two dimensional ones as La2-2£i;Sr2+2a;Mn207 the interplay of several de- 
grees of freedom, charge, spin, orbit and lattice displacements determine the physical properties of the system. The 
interactions between the different degrees of freedom cannot be reduced to perturbation theory as in other materials 
and the disorder produced by alloying to obtain the metallic state does not contribute to the clarification of theory 
nor to the interpretation of experiments. 

The recent observation of a insulator to metal transition (IMT) [1] in pure LaMnOa at 32 GPa and room temperature, 
well above the Neel temperature (145 K) and below the Jahn- Teller transition temperature (780K), opens the way to 
a study of the role of the orbital degrees of freedom on the electronic structure in a stoechiometric material. 

Two theoretical studies relevant to this matter have appeared in the literature almost simultaneously after the 
publication of the experimental result: both papers are based on local density approximation-t-Hubbard U (LDA-I-U) 
approximation. 

In the first one by Wei-go Yin et al.Q the electronic structure of the ground state of LaMnOa (antiferromagnetic 
A phase) is analysed to determine the relative importance of the electron-electron (e-e) against electron-lattice (e-1) 
Jahn- Teller interactions. It concludes that the e-1 interaction by itself is not sufhcient to stabilize the orbital ordered 
state and emphasizes the importance of the e-e interaction to facilitate the Jahn- Teller distortion. 

In the second paper by Yamasaki et al. Q LDA-f U and LDA+ dynamical mean field theories are used to analyze 
the metal insulator transition taking place at 32 GPa in the paramagnetic phase of the same material. The authors 
conclude that the transition at 32 GPa is caused by the orbital splitting of the Cg bands and that both e-e and e-1 
interactions are needed to explain the insulating character of the substance at lower pressures. 

In this paper we focus our attention on the orbital aspects of the insulator to metal transition. In order to include 
in an appropriate way the strong correlations induced by the dominant e-e interactions we follow Feiner and Oles in 
introducing auxiliary fields (Slave Bosons, SB) to the description of the low energy states. 

We start by defining a simplified Hamiltonian, two parameters: U on site Coulomb repulsion, and Ae a measure 
of the Jahn Teller splitting, to describe the tg states of the system and use SB to calculate the lowest energy varying 
the parameters to obtain a phase diagram. This procedure allow to identify the values of the parameters that are 
appropriate to describe the transition, as well as the thermodynamics and transport properties of LaMn03 in the 
paramagnetic phase. 
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II. METHODS 
A. Hamiltonian 



To describe the active electrons in LaMnOa we use the double-exchange model that contains the essential physics 
of manganiteslSj. The four 3d electrons in each Mn'^+ site are polarized in the same direction due to the large Hund 
coupling. Three of them occupy the t2g orbitals and are considered localized forming a spin S = 3/2, while the fourth 
occupying the eg state is itinerant. We treat the localized spin classically and since the fourth has to be parallel to 
the local spin we can consider the Cg electron as spinless with Hamiltonian 

H = Ht + Hu + HjT (1) 
where the first term, Ht, is the kinetic energy given by 

Ht-Yl ^f^lc,, (2) 

with a,/3 = 1,2 corresponding to the (possibly site dependent) orthogonal basis for the two Cg orbitals. The values 
of the hopping integrals t'^j^ depend both on the type of orbitals involved and on the direction between sites i,]. In 
LaMnOa, there is a staggered order in the x-y plane and the orbitals are stacked ferromagnetically along the z axis. 
The dominantly occupied orbitals alternating in the x-y plane are \x) = |3a;^ — r^) and \y) = |3j/^ — r'^'), which define 
then the low energy orbitals |1) in each of the sublattices, respectively A and B. The corresponding higher energy 
orbitals |2) are therefore respectively \y^ — z^) and \x^ — z^). The hopping parameters between these orbitals are the 
foUowings: 



^AB ^ . ^ 1/2 



-V3/2 
1/2 -V3/2\ 

y 



t being the hopping between \x) (jy)) orbitals along the x (y) direction. 

The magnetic order is introduced by modulating the hopping integrals by the factor exp(zAy) cos(6'y72),[l] with 
6ij being the angle between the t2g localized spins in the neighboring sites i,j and Aij a hopping phase. As we 
are interested in the paramagnetic phase that appears at room temperature T 3> T^v, we assume that the localized 
spins are completely random and consider a mean field approximation in which this factor is averaged giving a value 
(exp(iAy) cos{9ij/2)) = 2/3. The tf^f^ are the same than for a ferromagnetic phase with a factor 2/3. In the following 
we take this renormalized hopping as the reference t. 

The on-site Coulomb interaction between Cg electrons occuping both orbitals on the same site is given by 

Hu = uY^ Ui^na (4) 

i 

with riia = c'ia^ia the number operators. 

Finally, to model the effect of the Jahn- Teller (JT) deformation, we add a term that shifts the on-site energies of 
the Eg orbitals 1 . 2 in opposite directions 

HjT = Ae ^{ni2 - na) (5) 

i 

which corresponds to a JT splitting of 2Ae. 



B. Slave bosons method 



In order to treat the Hamiltonian Eq. we used the slave boson theory of Kotliar and RuckensteinQ adapted 
to our case of two orbitals instead of the two projection of spin. Therefore, we introduce new boson {ei,di,bia) and 
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pseudofermion (fia) operators. The boson numbers ejci, bl^bia and dldi represent the projectors onto the possible 
states |0i) (Oil, \ai) (a^l and \di) {di\ so that 



eje, + bl^bii + b\^ba + d\d, = 1 (6) 

and 

h + d'^ d- — r 

The original fermion operators are replaced by 

4i = {ble, + d\h,)fl 

4 = + (7) 

which correspond to a representation of the empty (|0i)), single occupied \2i)) and doubled occupied local 

states by 

|0.) = e\\vac) 
11.) = hyi\v^c) 
%) = byi\vac) 

\d^) = dlflfl \vac) (8) 

with \vac) correspondig to the vacuum state. 

The anticommutation rules for the original fermions are guaranteed provided the following constraints are satisfied 



+ - «-l-2 (9) 

which are implemented by means of the corresponding Lagrange multipliers {Ai, /^ii, ^^2}- To recover the correct 
result in the uncorrelated {U — 0) limit, a renormalization of the bosonic factor in Eq. ([7]) is necessary. In analogy 
with the spin case, the renormalized bosons factors take the form 



^ll'^i + dy 



'-il 



1 - e e 



y - blh,) (1 - did, - bib,,) 



z\, ^ . (10) 



b]nei + dlba 



1 - e\e, - b^by (1 - d\d, - b^ba 
With this slave boson representation, the Hamiltonian ([1]) reduces to 



A, 



(e\e, + blh, + b^ba + d\d, - l) (11) 



where the pseudo-fermions number operators are Uia = fJafia = cla'^ia- We have studied the solutions of this 
Hamiltonian in the mean-field approximation, in which we replace the boson operators by their averages obtained 
from the minimization of the band energy. In pure LaMnOa the number of conduction electrons is n = 1 and in the 
absence of charge ordering nn + 1112 = = 1 . In this case the band renormalization factor is independent of the type 
of orbitals involved q = zj^zjp. 
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III. RESULTS 



To characterize the solutions of Eq. pT|) for different vahies of the two parameters U and Ae, we considered both 
the magnitude of the pseudofermion band gap and the orbital polarization. In Fig. [1] we show the dependence of the 
band gap with the value of U. We can see that for Ae = there is a jump in the band gap at a critical value Uc, 
corresponding to a first order metal-insulating transition (in the parameter U). When we begin to increase the value 
of the JT splitting, Uc shifts to lower values and also the initial band gap in the insulating phase also diminishes. 
Beyond certain value of Ae (between 0.225^ and 0.250t), the metal-insulating transition is no more of first order, and 
the band gap opens smoothly. 



-A£/( = 0.0 
-Ae/« = 0.1 
-A£/t = 0.2 
■Ae/t = 0.3 
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FIG. 1: Dependence of the band gap as a function of U. 



Similar features can be observed in the orbital polarization shown in Fig. [21 where for small values of Ae there 
is a jump in the polarization at the same critical value Uc- From these results, we can see that an insulating phase 
with orbital polarization is possible even without a JT deformation. However, an orbital polarization will induce, 
due to the electron-lattice interaction, a JT deformation. We note also that for Ae = the orbital polarization is 
symmetric, i.e. the low energy orbital can be any combination of eg orbitals, and the actual one will depend on the 
electron-lattice interaction. 
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FIG. 2: Dependence of the orbital polarization as a function of U . 



In Fig. [3] we show the phase diagram U — Ae where the line corresponds to the value of Uc at the metal- insulating 
transition, and the gray tone represent the orbital polarization given by the occupancy of the low energy orbital n\. 
While in the insulating phase there is always a high orbital polarization, in the metallic phase the amount of orbital 
polarization depends on Ae: the higher the value of Ae, the higher the polarization. The critical value is close to 
Uc = 7t for the case without JT splitting, and descends to almost It for a splitting Ae = t. 

To better characterize the different phases we have calculated the pseudofermion density of states (DOS). In Fig.[3] 
we show this DOS for three different representative parameters in the phase diagram: 

% [U ^ 3.5t, Ae = O.Ot): metallic phase with very small orbital polarization which we have named orbital liquid 
(MOL). In this case, there is almost no orbital polarization (n\ — 0.502) and the renormalization factor is 
g = 0.917. 
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FIG. 3: Phase diagram [/-Ae. Solid line: metal-insulating transition. Gray tone: occupancy of the low energy orbital n\. 

• {U — 3.5i, Ae — QAt): metallic phase with orbital order (MOO). The orbital polarization is now higher 
{ni — 0.764) and the rcnormalization factor is closer to one [q = 0.941). 

• [U = 3.5i,Ae — 0.8t): insulating phase, also with orbital order (100). In this case there system is almost fully 
polarizad with m = 0.920 and also as a consequence q = 0.986 is also close to one. 
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FIG. 4: DOS for the different phases: (a) metallic phase with very small orbital polarization, (b) metallic phase with orbital 
polarization, and (c) insulating phase. 

To simulate the effect of pressure on the electronic properties, we take the experimental data on the IM transition 
and model the dependence of the parameters with the pressure. All magnitudes are given in terms of the effective 
hopping at P = 0, which is of the order of to ~ 0.4 eVQl. The JT splitting is taken to vary linearly from a value of 
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Ae = 1.85io[l] at P = to Ae = at the experimental JT suppression pressure P = 18 GPa. The e-e interaction 
constant is considered constant. FinaUy, the effective hopping has a dependence of the form[2|] 

where d(P) is the mean Mn-0 distance as a function of the pressure, taken from the experimental data, and do = d(0). 



"1 

0.5 0.6 0.7 0.8 0.9 1 




P(GPa) 

FIG. 5: Phase diagram U-P. Solid line: metal-insulating transition. Dot line: U = ll.Sto- Gray tone: occupancy of the low 
energy orbital ni. 

From this dependence of the parameters with the presure, in Fig. O we show the phase diagram U — P where as 
before the solid line corresponds to the value of Uc at the metal-insulating transition, gray tone represent the orbital 
polarization given by the occupancy of the low energy orbital ni . Note that the values of U in the vertical axis are 
in units of to, and we have added a dot line showing that for the estimated value of U = 11.5to we obtain a insulator 
to metal transition pressure close to the experimental one of P ~ 32 GPa. 

IV. CONCLUSIONS 

We use a minimum parameters model Hamiltonian to study the evolution of orbital polarization in LaMnOa. In 
order to include appropriately the effects of correlation, we resort to the Slave Bosons technique which was previously 
used by Feiner and Oles in the context of manganites. 

We calculate the electronic structure and from it, the electronic energy to obtain a phase diagram in terms of two 
independent parameters U/t an Ae/t. The results can be translated to the effect of pressure on the material by 
modelling the variation of the parameters with volume and connecting to pressure though the compressibility. 

The same Hamiltonian could be used to represent the orbital state of other compounds where trivalent Rare Earths 
substitute partially or totally La, again though modelling of the variation of the hopping , U or Ae as an effect of 
substitution and pressure, as for example in equation 1121 

Our results confirm and make transparent the conclusion reached in previous ab-initio calculations: The inclusion 
of the Coulomb energy is necessary and constitute an important factor enhancing the orbital polarization in these 
compounds. 

From the density of states it is possible to calculate the gap in the insulating phases and the number of carriers as 
a function of temperature and pressure, in order to compare with the results of Loa et al.|l| 

A natural continuation of these results is the evaluation of the elastic energy involved in the Jahn- Teller distortion 
that will appear at the polarized phases, the results that will be published elsewhere. 
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